Relics of cosmological quark-hadron phase transition 
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We propose that the amplified density fluctuations by the vanishing sound velocity effect during 
the cosmological quark-hadron phase transition lead to quark-gluon plasma lumps decoupled from 
the expansion of the universe, which may evolve to quark nuggets (QNs). Assuming power-law 
spectrum of density fluctuations, we investigate the parameter ranges for the QNs to play the role of 
baryonic dark matter and give inhomogeneities which could affect big-bang nucleosynthesis within 
the observational bounds of CMBR anisotropy. The QNs can give the strongest constraint ever 
found on the spectral index. 



PACS numbers: 98.80.Cq, 12.38. Mh, 95.35.+d 

As the temperature of the universe cools down to 
the critical temperature Tq ~ 100 — 200MeV, so called 
quark-hadron phase transition occurs in such a way that 
the high temperature quark gluon plasma (QGP) phase 
transforms to the low temperature hadronic phase. If the 
phase transition is first order, it is described by the nu- 
cleation of hadron bubbles and their growth. As hadron 
bubbles occupy more than half of the whole space, the 
picture becomes that of QGP bubbles shrinking in the 
sea of the hadronic bubbles. Witten pointed out that the 
shrinking QGP bubbles in the transition may evolve to 
quark nuggets (QNs) which can play the role of baryonic 
dark matter [Q. Baryon concentration in the shrinking 
QGP bubbles also give inhomogeneities which could af- 
fect big-bang nucleosynthesis (BBN) Even though 
recent lattice calculations show that the quark-hadron 
phase transition is weakly first-order 0, the mean bub- 
ble separation lb ~ 1cm based on the lattice calculations 
PJ is too small for the transition to achieve the effects 
mentioned above. 

Recently, it was shown that the growth of the subhori- 
zon scale fluctuations in the quark-hadron transition is 
amplified because the sound velocity during the transi- 
tion vanishes ||[|. Amplified subhorizon scale fluctua- 
tions give gravitational potential wells which will trap 
cold dark matters (CDMs) like axions and primordial 
black holes (PBHs) Q . Resulting CDM lumps have size 
of about 1 — 10AU today. Also the vanishing sound veloc- 
ity effect enables, though it needs fine-tuned initial con- 
ditions 0, the formation of PBHs of 1M Q ~ 10 33 g which 
can be related to MACHOs [p[. The amplified fluctua- 
tions may produce nonlinear lumps which are decoupled 
from the expansion of the background universe. If the 
decoupled lumps were hadronic, however, they would be 
washed out by neutrino damping || before the big-bang 
nucleosynthesis era [||J^] , so could not give any significant 
effect on dark matter problem as well as the inhomoge- 
neous big-bang nucleosynthesis (IBBN) model. In this 
letter, we propose that the fluctuations amplified by the 
vanishing sound velocity effect may produce nonlinear 
QGP lumps which will survive till today in the form of 



QNs. With a power-law assumption on the density fluc- 
tuation spectrum, we find possible parameter ranges for 
which the QNs form to play the role of dark matter and 
to produce baryon number inhomogeneities for the IBBN 
within the observational bounds of CMBR anisotropy. 

The sound velocity vanishes during the phase transi- 
tion as follows. The QGP and the hadronic phase are in 
phase equilibrium with constant temperature Tq shortly 
after reheating. In the cosmological situation, chemical 
potential of each fluid is negligible. Thus the pressure 
and the energy density depend only on the temperature 
of the system. Consider the system with a scale smaller 
than the horizon size at the transition d^j ~ 10km but 
larger than the mean bubble separation lb, where the de- 
tails of bubble dynamics are not important. Before and 
after the transition, the system is filled with QGP and 
hadronic radiation fluid, respectively, so the sound veloc- 
ity v s = -\/l/3. As the phase transition proceeds, while 
the pressure remains constant because the two phases 
are in pressure equilibrium, the energy density of the 

system gradually decreases. Therefore, the sound ve- 

l li 

locity v s — (dp/dp) s = during the phase transition 
Since the sound velocity vanishes, the pre-existing 
density fluctuations grow without any pressure gradient 
and restoring force and the amplitudes of the fluctua- 
tions are amplified. Superhorizon scale fluctuations are 
not affected by the sound velocity transition. 

The subhorizon fluctuations which entered the horizon 
long before the quark-hadron transition were damped by 
the neutrino diffusion. Only the fluctuations with the 
comoving wavenumber k > 10 4 fc# can retain their hori- 
zon crossing amplitude Sh , where ku is the comoving 
wavenumber of the fluctuation which enters the horizon 
at the transition. During the whole time of the transition, 
the sound velocity transition amplify the amplitudes of 
the subhorizon fluctuations with and k upto 

S a H mP =(£")' Sh for k H^k<ltfk H (1) 
where £ = 1 for the bag model and £ = 3/4 for lattice cal- 
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culation The exponent ( is obtained for k^> k B but 
we simply regard it valid for k > k B - The fluctuations 
with 5 H rnp > 1, arc nonlincarizcd by the sound velocity 
transition and are decoupled from the expansion of the 
background universe. So it needs 8 B ^ 5 C H = (k B /k) 1 ^ for 
the nonlinearization. Since this condition is far weaker 
than 5 C H ~ 1 in ordinary case without the sound veloc- 
ity transition, the nonlinear objects with mass below the 
horizon mass at the transition w 10 33 g (even PBHs 
with 10 33 g ||) are formed much more copiously during 
the quark hadron phase transition. 

What will happen as the density fluctuations grow and 
become nonlinear? We propose that some nonlincarizcd 
fluctuations lead to the formation of QNs. The non- 
linearization of fluctuations with scale A in the linear 
analysis can be interpreted in another sense as the over- 
dense regions of size A evolving into gravitationally bound 
lumps decoupled from the expansion of the universe. The 
temperature of the decoupled lumps will not decrease any 
more by the expansion of the universe. Instead, the in- 
creased pressure and energy density of the lumps by the 
nonlinearization (or the gravitational binding) will in- 
duce the transition to the QGP phase in the hadronic 
phase regions. If the overdense regions were composed 
of a single phase, the temperature of the overdense re- 
gion, in general, would increase as the energy density 
increases. For a region with mixed phases of QGP and 
hadronic phases, however, the temperature of the region 
does not increase until the hadronic phase is converted 
entirely to the QGP phase. The phase of the nonlin- 
earized lump is determined by comparing two important 
but uncertain time scales. One is t non ii the time needed 
for the overdense region to become nonlinear. The other 
is tback, the time for the mixed phase to return to QGP 
phase. What is needed for the decoupled lumps to return 
to the QGP phase is t non i ~ t back . If t back > t nonh the 
decoupled lumps are made of almost hadronic phase so 
that they will be washed out by the neutrino damping. 
Also, if tback < tnonU the lump is still in linear regime 
when it is entirely converted back to the QGP phase and 
it experiences ordinary phase transition to the hadronic 
phase due to the cooling by the expansion of the uni- 
verse before it becomes nonlinear. Therefore, it needs 
tnoni ~ tb ac k for nonlinear QGP lump formation. This 
can be achieved for the overdense regions with 5h ~ 5%- 
For such regions t non i ~ A£q that is, roughly the time 
needed for the phase transition between two phases. So 
we can conceive that the overdense regions with 8 B ~ 5 C H 
evolve to nonlinear QGP lumps. Though the overdense 
regions with 8 B 5 C H which have t non \ AtQ may also 
satisfy the condition t non i ~ tback, the number of such 
regions is negligible as can be seen in the probability dis- 
tribution of overdense regions. 

The mass and size of the QGP lump are specified by 
the baryon number in the lump N B oc k~ 3 as Mq = 

{ J^b ) M Q where N B - H = 10 49 (T Q /100McV)- 3 is the 

total baryon number contained in the horizon at the tran- 



sition. For k « 10 4 fcff, the minimum scale relevant to 
the sound velocity transition, N B ~ 10 37 . The QGP 
lump can evolve to a QN if the quarks in the lump 
become degenerate. It occurs if the baryon density to 
entropy density ratio, n B /s ~ 1. For comparison, the 
background universe has (ng/s)^ ~ 3.8 x 10~ 9 (tt B h 2 ) 
where VL B is the density fraction of the baryonic mat- 
ter at present and h is the Hubble parameter in units 
of 100km sec" 1 Mpc -1 . The resulting QN mass M QN 
is determined by how and when the lump get enormous 
increase in n B / s and transform to a QN. Here, we do 
not attempt to answer these questions but will give some 
qualitative descriptions. For this we introduce a param- 
eter k, the QN to the QGP lump volume ratio. If one 
takes the quark matter energy density ~ the QGP en- 
ergy density at the transition, then Mqn = kMq. From 
a chromoelectric flux-tube model (CFTM), it was found 
that the baryon penetration at the surface of the QGP 
lump is very low and the baryon number density in the 
lump increases as the surface hadronizes [ |io| |. With this 
view, QNs may form by increased ns during the sur- 
face hadronization JlO[ ] . If there were no efficient ways 
to get rid of entropy and the QN formed solely by the 
surface hadronization (only by the increase of n B ), then 
K ~ {n B /s)bg and ns should be increased to about the 
baryon number density of a nucleon ~ 0.3fm~ 3 . How- 
ever, for the QGP lump with N B < 10 34 (n B /s) fcg 1 ~ 

10 42 (fis/i 2 ) -1 ) the lump becomes smaller than the neu- 
trino mean free path l v ~ 10(T/100MeV) _5 cm before 
reaching k ~ (n B /s)b g - Then the entropy in the lump 
will be reduced rapidly by neutrinos. If this occurs faster 
than the surface hadronization, then k can be approx- 
imated as k ~ (ATg/lO 34 ) . In fact, the entropy loss 
by neutrinos also occur at distances < l v from the sur- 
face and it needs more detailed calculations to find the 
condition for QN formation. However, we simply assume 
k ~ (Nb /lO 34 )^ 1 as the upper limit for the QN forma- 
tion. Resulting QNs have masses < 10 18 g. In the models 
generating large entropy before the quark-hadron transi- 
tion jn| , the QNs can be formed easily with k ~ 1. This 
situation is, however, not relevant to our work because 
there can be no sound velocity transition. 

Once the QNs form, they are unstable to surface evap- 
oration [jl2] and boiling Only the QNs with the 
baryon number larger than the critical baryon number 
N B can survive till today. With small baryon penetra- 
bility in the CFTM, N B is lowered and two results of 
N C B = 10 39 , 10 44 are known @Q. However, as claimed 
by the authors, their works overestimated the baryon 
evaporation. Especially, they assumed flavor equilibrium 
in their calculations. The flavor non-equilibrium will re- 
duce N B . So, we take N B as a free parameter bearing in 
mind the above values. The baryon number of the QNs 
needed to affect the IBBN is smaller than N B , but the 
difference will not be large considering the rapid evapo- 
ration rate below N B HO. 

Now, we have found the conditions for the QN forma- 
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tion, Sjj ^ Sjj and Nb ft N B . To estimate the number 
of the QNs, we just count how many overdense regions 
satisfying the conditions exist. We would like to empha- 
size that the number density of the QNs depends much 
on the details of density fluctuation. Though the details 
of the phase transition may rule out the QN formation 
from the shrinking QGP bubbles, whether the QNs from 
the sound velocity transition can be produced enough 
is another problem. One thing to note further is that 
the QN formation proposed here occurs even without the 
sound velocity transition but with much severe condition 
Sh }t 1- The sound velocity transition relax the condition 
and enables more copious QN formation. 

We assume a simple power-law spectrum of density 
fluctuation with \Sk\ 2 oc k n , where 8k is the Fourier trans- 
form of <5(x) = (p(x) — pb)/pb and the spectral index n 
is a constant. The initial power spectrum of the fluc- 
tuation amplitude is defined by the rms amplitude for 
a given logarithmic interval in k, 5 2 (k) = fc 3 |<5fc| 2 /27r 2 . 
From the linear analyses, the rms amplitude 5~i(k) at ti, 
the time when the fluctuations develop, is related to the 
rms horizon crossing amplitude 5n{k) as follows 



S H (k) 



k 
km 



Si(k) 



(n-l)/2 



S H (k ) (2) 



where km is the wave number of the horizon scale at ti 
and the subscript '0' denotes the values at present. From 
the COBE measurement, Sii(ko) ~ 10~ 5 and n — 1.2 ± 
0.3 fl6f| . The number density of overdense regions with 
Sh ^ 8 C H can be found by the Press-Schechter method 
|l7f which is widely used to estimate the number density 
of structures like galaxies. The difference here is in the 
scale dependence of the critical amplitude S H oc fc~^. 
The initial mass spectrum of the QGP lump in the range 
(M H , Mh+cLMr) where M H oc Mq /3 is the horizon mass 
of scale with k 
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where a = n + 2C~l and a H = a Ho (M H / ' M Ho ) ( - 1 - n '>/ i is 
the filtered amplitude of <5f/(/c) and ohq = 9x 10~ 5 [ 18l . 
The cutoff scale M B <c M H . Due to the scale depen- 
dence of S H , the mass spectrum has exponential cutoff 
even with n ~ 1 Harrison-Zel'dovich spectrum. For com- 
parison, n(M H ) oc Aftf ( ™~ 11)/4 expJ-M ff (n - 1)/2 ) with- 
out the sound velocity transition pjj ], so if n = 1, the 
mass spectrum could have very broad mass ranges with 
n(M H ) oc M H ~ 5/2 . With the exponential cutoff, the QN 
with only Nb — N B form significantly. So we can as- 
sume 5-function type mass spectrum and approximate 
the initial number fraction of the QNs with N B as 



P i (Mf I )=a H (M H )exp 



'H 



V2a H {M- H ) 



(4) 



which is used for calculating the number fraction of PBHs 
formed by the density fluctuations with blue-shifted spec- 
trum. 

Since the QNs can be regarded as pressureless dust, 
the density fraction at present is 
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&QN(to) = &>QN(teq) = 



Peq 



M, 



QN 



a(U) 



A (5) 



where the subscript 'eg' represents the values at matter- 
radiation equal time. M c is the mass contained in the 
overdense region with k c at ti. The density fraction can 
be arranged into 



^QN(t ) = K 



eg 



Pi 



(0) 



Fig. 1 shows the relation between k and N B to satisfy 
Oqat = 1 for possible range of the spectral index, 0.9 < 
n < 1.5. It seems impossible for QNs to be formed solely 
by the surface hadronization (k ~ (tib / 's)bg)- The QGP 
lumps become smaller than l v when k = (iV|/10 34 ) _1 
(the bold line in Fig. 1). So, in the left side of the bold 
line, the QGP lumps can evolve to QNs (ns/s ~ 1) by 
further losing their entropy by neutrinos. The QN for- 
mation is relevant for only iV^ < 10 42 . For n « 1, it 
needs N B < 10 38 for £ = 1, which is smaller than the 
CFTM results With ( = 0.75, the QN formation 

is possible if n > 1.1. Assuming k = (Ng/10 34 )^ 1 , the 
upper limits on n are found from the condition Qqn < 1 • 
If N B < 10 40 , the QN can give strongest constraints ever 
found, 0.9 < n upper < 1.2 (the bold lines in Fig. 2). 
PBHs can give at best n upper > 1.23 even weaker as the 
reheating temperature increases |l8| , p0|| . Note that the 
upper limits on n do not depend on ti or the reheat- 
ing temperature in inflationary models. Demanding less 
Qqn, the parameter ranges become broader and the con- 
straints on n becomes stronger. 

Since QNs have positive electric surface potential of or- 
der ~ MeV, they absorb only neutrons. This reduces the 
neutron to proton ratio, so lowering the helium produc- 
tion. Not to violate the observations, the QNs should 
have size ^ 10 _6 cm assuming £Iqn = 1 and perfect 
baryon penetrability (even smaller with smaller penetra- 
bility) HI. The corresponding N% ~ 10 13 , so the QNs 
considered here easily satisfy the condition. More signifi- 
cant effects on the IBBN can be induced by the evapora- 
tion of the QNs. The minimal requisition for the IBBN is 
that the mean separation between QNs should be larger 
than proton diffusion length when the BBN starts. It 
requires the mean separation of QNs at the transition, 
Iqn ^ 3m ruling out the QNs from the shrinking 
QGP bubbles (4 w 1cm). We find 
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The upper values of n for the IBBN (Iq N ss 3m) are 
shown in Fig. 2. It can be seen that the the mini- 
mal condition is easily satisfied. The closure condition 
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Qqn < 1 (the bold lines in Fig. 2) corresponds to de- 
manding Iqn ^ 50m. It goes to Iqn > 100m demanding 
£Iqn < 0.1 not shown in Fig. 2. Interestingly, demand- 
ing 0.1 < £Iqn ^ 1, Iqn lies in the ranges for the IBBN 
to be effective 
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So, the QNs can contribute the 
density of the universe in the form of baryonic dark mat- 
ter as well as the IBBN. The enhanced heavy element 
formation can be the signature for the QNs | f23| . 

In summary, we propose that QNs can be formed by 
the sound velocity transition during the quark-hadron 
phase transition. The difficulties of the QNs formed 
by the shrinking QGP bubbles can be avoided by the 
QNs from the sound velocity transition because their for- 
mation depends much on the details of density fluctua- 
tions. So, they can dominate the density of the universe 
Oqjv ~ 1 and affect the standard BBN. Also, the QN 
formation can strongly constrain the spectral index. Our 
analyses are so far rather qualitative and include unde- 
termined parameters such as k and Ng . It needs further 
systematic analyses to get more quantitative results. Af- 
ter all, the sound velocity transition effect can make the 
QN to play an important role during the evolution of the 
universe. 
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FIG. 1. The parameter ranges for « and Ng within the 
observational bounds of n — 1.2 ± 0.3. for Tq = 150MeV 
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FIG. 2. The relation between the spectral index and the 
inhomogeneity scale Iqn for Tq — 150MeV. The bold 
lines are the upper limits on n found from Qqjv = 1 with 
k = (W 34 /N B ) 
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